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A detailed derivation of the master equation of the cavity radiation of a coherently prepared 
y-shaped four-level correlated emission laser is presented. The outline of the procedures that can 
be employed in analytically solving the stochastic differential equations and the rate equations of 
various correlations are also provided. It is shown that coherently preparing the atoms in the 
upper two energy levels and the lower, initially, can lead to a genuine continuous variable tripartite 
entanglement. Moreover, preparing the atoms in the coherent superposition, other than the possible 
maximum or minimum, of the upper two energy levels, leaving the lower unpopulated, may lead 
to a similar observation. With the possibility of the atom at the intermediate energy level to take 
three different transition roots guided by the induced coherence, this system in general is found 
to encompass versatile options for practical utilization. In particular, coupling at least one of the 
dipole forbidden transitions by an external radiation is expected to enhance the degree of detectable 
entanglement. 

PACS numbers: 42.50.Ar, 42.50.Gy, 03.65.Ud 



I. INTRODUCTION 

In recent years, due to the relative simplicity and high 
efficiency in the generation, manipulation, and detection 
of optical continuous variable (CV) states [JJ, the cor- 
responding entangled states have been successfully im- 
plemented in unconditional quantum teleportation [3] , 
quantum dense coding [4], quantum error correction [5], 
and universal quantum computation [DJ among others. 
With the progress in the CV entanglement research, the 
generation of more than a bipartite entanglement has 
attracted much attention. A truly A-partite entangled 
state generated by a single-mode squeezed state and lin- 
ear optics [7J along with the generation of CV tripartite 
entanglement using cascaded nonlinear interaction in an 
optical cavity without linear optics have been theoret- 
ically investigated [5] and also experimentally realized 
[DJ. Nevertheless, the structure of the entanglement 
for the three-mode system is a bit more subtle than that 
for a bipartite case, wherein, different classes of entangle- 
ment are defined based on how the density matrix may 
be partitioned [TO] . The classifications range from fully 
inseparable, which means that the density matrix is not 
separable for any grouping of the modes (genuine tripar- 
tite entanglement), to fully separable, where the three 
modes are not entangled in any way. Despite the chal- 
lenge, fortunately, there is a large number of works that 
are aimed at devising ways of detecting a tripartite en- 
tanglement [TUrfi!] . 

On the other hand, there has been enormous effort in 
studying quantum optical systems that are capable of 
generating tripartite entangled light [5J [TOJ [H] that in- 
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eludes the nondegenerate parametric oscillation and six- 
wave mixing in the nonlinear medium. As an alternative, 
the idea of generating a tripartite entanglement from var- 
ious schemes of four-level atomic systems via coherent 
superposition has been under study recently [16f|22| . En- 
tangled photons from these systems are expected to have 
a potential applications in quantum memory |23j and 
long distance quantum-communication |24j , since the low 
frequency and narrow linewidth of the light can enhance 
efficient coupling between photons and atomic memories 
in a quantum network [23 26 . Some of the possible 
schemes include but not limited to A- type [16], V-type 
[17] , cascade-type PH HH], and Y-type [2D]. In these 
works, the coherent superposition is induced by excit- 
ing an atom in a lower energy level to the upper (using 
an external pumping mechanism) from where the atom 
undergoes direct spontaneous emission that leads to a 
genuine tripartite entanglement [TBI [2J • Quite recently, 
Shi et al. [20] extensively discussed the way of generating 
a tripartite entanglement applying a mechanism of six- 
wave mixing. The three external radiations they applied 
were believed to be responsible for creating the required 
nonclassical correlations. 

However, in this contribution, following a similar kind 
of reasoning, it is proposed that the V-shaped four-level 
atomic scheme can be a reliable source of strong entan- 
gled light if the initial preparation of the atoms is as- 
sumed to be accountable for inducing coherent superpo- 
sition. Moreover, rather than placing the atoms in the 
cavity throughout the operation and expose them to ther- 
mal fluctuations, as usually the case, it is assumed that 
they are injected into the cavity at a constant rate. For 
the sake of convenience, the amplification of the light 
due to reflection between the walls of the coupler mir- 
ror when photons with different frequency emitted in the 
forked four-level cascade transition are correlated by the 



2 



coherence induced via initial preparation is dubbed as 
a nondegenerate ^-shaped four-level correlated emission 
laser. Since a large number of atoms, in principle, can 
be injected into the cavity over a longer period of time 
and the direct spontaneous emission process can also be 
efficient, if the atoms are properly prepared initially, it 
would be reasonable expecting this system as a reliable 
source of bright light. 

In order to study the dynamics of the entanglement 
employing the existing criteria, it is found necessary and 
appropriate establishing the mathematical framework for 
solving the involved differential equations beforehand. To 
this effect, based on the involved structure of the atomic 
levels, the master equation is derived following the out- 
line presented elsewhere for the corresponding two-mode 
case [27J HFJ ■ Due to the coherent superposition induced 
via initial preparation and cascading process, a signifi- 
cant correlation in the generated three modes is observed 
in the calculated master equation. Moreover, with the in- 
tention of paving the way for in depth analysis, the pro- 
cedure of solving the emerging coupled differential equa- 
tions is outlined. It is found that the required solutions 
can be explicitly written down once the corresponding 
3x3 matrix constructed from the prefactors in the mas- 
ter equation is diagonalized and its cigcnmatrix is con- 
structed, which in principle is a surmountable task al- 
though the rigor may be somewhat lengthy. 



II. DESCRIPTION OF THE MODEL 

In earlier studies on the three-level cascade laser, it 
was observed that the spontaneous transition during the 
cascading process induces a coherent superposition that 
leads to enhanced quantum features including a bipar- 
tite CV entanglement [29l [30] . It can be asserted that 
initially preparing atoms in a certain coherent superpo- 
sition and then allowing them to follow realistic spon- 
taneous transition roots can yield a strongly entangled 
light [31T433J . Taking this as a motivation, the y-shaped 
four-level atomic system initially prepared in a coherent 
superpostion of the energy levels between which a direct 
electric dipole transition is forbidden would be consid- 
ered. For the sake of convenience, the lower energy level 
is denoted by |0), the intermediate energy level by |1), 
and the upper two energy levels by |2) and |3) (Please 
note that the schematic representation of the involved 
atomic energy levels is provided in Fig. [I]). In order to 
expedite the cascading process, it is assumed that the 
parity of the energy levels |0), |2), and |3) is the same 
whereas that of |1) is different. This entails that direct 
spontaneous transitions between energy levels |2) •<->• |3), 
1 2) «-» |0), and |0) f-> |3) are electric dipole forbidden, but 
due to the parity difference, the transitions between |I) 
O |0), |2) O |1), and |1) <H> |3) are allowed. It is worth 
noting that if required the dipole forbidden transitions 
can be induced by an external pumping mechanism in a 
similar manner as in the three- level case |34H37j . 



While the atom undergoes a direct spontaneous tran- 
sition from energy level |3) to |1), suppose it emits a 
photon represented by annihilation operator S3. In prin- 
ciple, it can still undergo a direct spontaneous emission 
and go over to the lower energy level |0); in the process 
emits a photon described by d\. In the cascading tran- 
sition from energy level |3) to |0) via |1), a correlation 
between the two emitted photons (03 and di) can readily 
be established. In a similar manner, it is not difficult 
to realize that there could be a correlation between the 
photons emitted (S2 and Si) when the atom undergoes a 
spontaneous transition from energy level |2) to |0) via 
These two processes, which are not entirely independent, 
are expected to initiate nonclassical correlations between 
the photons emitted while the atom cascades from the 
upper two energy levels to the lower via different forked 
roots. Nevertheless, to establish a genuine tripartite en- 
tanglement as prescribed by the von Loock and Furusawa 
criteria a correlation between the photons emitted 
from the upper two energy levels (S2 and S3) is very cru- 
cial. In order to initiate this important correlation, it is 
worth noting that in case the atom is initially prepared in 
a coherent superposition of the upper two energy levels 
and the lower, it does not arbitrarily undergo the afore- 
mentioned spontaneous transitions due to the resulting 
population sharing. 

If there is a triply resonant light in the cavity, the 
atom with energy level |1) has three distinct alternatives 
to take except for the spontaneous decay to any other 
energy level that is not involved in the present consid- 
eration. The first and the most natural one is to con- 
tinue with the direct spontaneous emission and then goes 
over to the lower energy level, or it can absorb a photon 
(S2) and excited to the upper energy level |2), or it can 
absorb a photon (S3) and excited to the upper energy 
level |3). In this description, if one manages to send 
a large number of initially prepared atoms through the 
cavity, the absorption-emission mechanism painstakingly 
follows different roots, which leads to additional correla- 
tion between emitted photons. In this scenario, as long 
as there is a coherent superposition between the energy 
levels 1 2) and |3) initially, a meaningful correlation in a 
subsequent emission of photons denoted by S2 and S3 is 
expected. It is, hence, envisaged that this process can 
forge the required nonclassical correlation between the 
photons emitted from the upper two energy levels. 

The more general explanation and possible setup for 
practical utilization of similar schemes were provided ear- 
lier in Refs. [351 131] • However, this contribution has 
one essential difference in which the initial preparation is 
assumed to be a prominent source of the atomic coher- 
ent superposition as opposed to the external pumping 
mechanism. Although the initial preparation and injec- 
tion process lead to some technical difficulties in practical 
utilization of the potential of this system, it is envisaged 
that the challenge would be less severe when compared to 
an external pumping mechanism that naturally initiates 
thermal fluctuations and atomic broadening. 
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from imperfect preparation is not taken into considera- 
tion. 

The atoms prepared in this manner are assumed to be 
injected into a triply resonant cavity at a constant rate 
r a and removed after sometime T, which is long enough 
for the atoms to spontaneously decay to energy levels 
that do not involve in the process. The density operator 
for the cavity radiation plus a single atom injected into 
the cavity at time tj can be denoted by PAR,(t, tj), where 
t — T < tj < t. The density operator for all the atoms 
in the cavity plus the cavity radiation at time t can be 
expressed as 



PAn{t) 
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PAFt(t,tj)Atj 



(4) 



FIG. 1: Schematic representation of the nondegenerate Y- 
shaped four-level atom. It is assumed that the atom under- 
goes direct spontaneous transitions from energy level |3) to |1) 
in a process emits a photon describable by annihilation oper- 
ator as, when it goes from |2) to |1) ai, and when from |1) to 
|0) oi. The three transitions are presumed to be resonantly 
coupled to the cavity radiation. 



where r a Atj represents the number of atoms injected into 
the cavity in a time interval of Atj. Assuming that the 
atoms are continuously injected into the cavity and tak- 
ing the limit that Atj — > 0, the summation over j can be 
converted into integration with respect to t' , 



PAn{t) 



p AR (t,t')dt'. 



(5) 



t-T 



III. MASTER EQUATION 

The interaction of a nondegenerate l^-shaped four-level 
atom with a triply resonant cavity radiation can be de- 
scribed in the rotating-wave approximation and the in- 
teraction picture by the Hamiltonian of the form 



H I =ig[a 3 \3)(l\-\l)(3\al 
+ a 1 |l)(0|-|0)(l|a t 1 ], 



a 2 |2)(l|-|l)(2|a T 2 



(1) 



where g is a coupling constant chosen to be the same for 
all transitions for convenience and a^s are the annihila- 
tion operators that represent the three cavity modes. 

Assuming that the atoms are initially prepared in the 
coherent superposition of the atomic energy levels except 
the intermediate, the pertinent atomic state can be taken 
as 

|*a(0)) = C 3 (0)|3) + C 2 (0)|2) + C (0)|0>, (2) 

where Ci (0)'s are the probability amplitudes for the atom 
to be initially in the i's energy level and the correspond- 
ing density operator takes the form 



P$|3><3|- 
+ /4° 3 ) |0)(3|- 
+ 4°o ) |0}(0|, 



-P$|3><2| 
P2° 2 ) |2>(2|- 



-P$|2><3| 



■Po°2 ) |0)(2| 



(3) 



where p^P's are the initial populations and pjj°L-y)' s are 
the coherences between the atomic energy levels. It is 
worth noting that the intermediate energy level is ini- 
tially unpopulated and the phase fluctuation resulting 



Replacing the summation over randomly injected atoms 
to integration in a similar manner has been done fre- 
quently [571I1D1I1I]. 

It is also a well established fact that the density oper- 
ator evolves in time according to 



d_ 



pA R {t,t') = -i[H, p A R(t,t')]. 



(6) 



It is not hard to observe that t' can be switched in such a 
way that pAn(t, t) represents the density operator for an 
atom plus the cavity radiation at a time when the atom is 
injected into the cavity, whereas PAn(ti t — T) represents 
the density operator when the atom is removed from the 
cavity. Since the atomic and radiation variables are not 
correlated at the instant the atoms arc injected into or 
removed from the cavity, it is possible to propose that 



PAR(t,t) = p A (0)p(t), 



p AR (t,t-T)=p A (t-T)p(t)> 



(7) 



(8) 



where p A (0) = pA{t). 

Hence, in view of Eqs. ^ and Q, integration of Eq. 
§5§ results 



dt 



PAR(t) = r a [p A (0)-p A (t-T)]p(t)-i[H, PA R {t)]. (9) 



Now taking the trace over the atomic variables using the 
fact that Tr A {p A (0)) = Tr A (pA(t — T)) = 1, leads to 



dp{t) 
dt 



-iTr A [H, pAn(t)}. 



(10) 
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Upon employing Eqs. and (10 1, the time evolution 
of the reduced density operator for radiation turns out 
to be 



dm 

dt 



■ s[P3l4 - 4^31 + «3Pl3 - Pl3«3 

P2ia\ - a\p 2 \ + d 2 p 12 - Pu^s] 



Pwa\ 



(11) 



in which = (i\pAFi\j) with i, j = 0, 1, 2, 3. 

On the other hand, on the basis of Eq. ffify, one can 
readily write 



d_ 
dt 



Pij(t) = r a (i\p A (0)\j)p - r a (i\p A (t - T)\j)p 

-i{i\[H, p A R{t)]\j)-lPij, (12) 



where the last term is introduced in order to account for 
the atomic decay process. Basically, 7 is the decay rate 
associated with every atomic transition including the rate 
of dephasing. In practical situation, assuming all decay 
rates as equal may not be reasonable as recently discussed 
elsewhere [251152] . 

Assuming the atoms to be removed from the cavity af- 
ter they have decayed to energy levels that do not involve 
in the lasing process implies that (i\pA(t — T)\j) = 0. 
Hence, on account of Eqs. ([I]), ([3]), and (12), one gets 

d „ , . 



IPij + r a p(t) [p 33 S i3 8 3j + P32 <*i3* 



(0); 



(0)r r . (0) jr X , (0); j (0) r r 

+ P23 "i2"3j + P30 0i3O j + p Q3 Oio03j + p 22 Oi 2 2 j 
+ PwfaSoj + Po2^i0^2j + Poo Sio5oj] 

+ g[a-3Pi]Si3 - alp 3j Sii + a 2 pijS a - alfajSn 
+ axpojSu - d\p lj S 0l - pi 3 a 3 8^ + paa\8 3i 



pi2a 2 Sij + Pii&lS 2j - pnaiSoj + Pio&lSxji, 



from which follows 



P3z(t) = r a pf z ] p{t) + g[a 3 pu + P3ia j 3 ] - 7P33, (14) 



dt 



d „ , . 



r a PyiP{t) + g[a 3 p 12 + p 31 d 2 } - 7p 32 , (15) 



d 
dt 



p2o(t) = r a p$o p{t) + g[a 2 p 10 ~ p 2 iOi] - JP20, (20) 



■^Pn(t) = gfiiPoi - /5i3«3 - P\2d 2 + p w a\ 
- a\p 3 i - a\p 2 {\ - 7Pii) 



(21) 



^PioM = g[aiPoo - Piih ajpso - 4p2o] - 7Pio, 



(22) 



Poo(*) = r a PooP{t) - g[dip w + p 01 a[] - -/p 00 . (23) 



In the good cavity limit (7 3> K, where k is the cav- 
ity damping constant), the cavity mode variables change 
slowly compared with the atomic variables. It is, hence, 
expected that the atomic variables will reach steady state 
in a relatively short time. The time derivative of such 
variables can be set to zero, while keeping the remain- 
ing atomic and cavity mode variables at time t. This 
procedure is referred to as the adiabatic approximation 
scheme. Confining to linear analysis, which amounts to 
dropping the terms containing g in Eqs. (14 1, (15), (17), 



(18), (20), and (23), and applying the adiabatic approx- 



imation scheme, one finds 



Pij 



(0) 
r aPij 



Pit), 



hi = 0) 



(24) 



(25) 



(13) with p = p(t) and i,j — 0,2,3. Eq. (25) indicates the 



procedure in which the initially unpopulated energy level 
|1)) is adiabatically eliminated. 



At this juncture, making use of Eqs. (16), (19), (21), 



(22), (24), (251, and applying the adiabatic approxima- 



tion scheme once again, it is possible to verify that 



P31 



gr a p 
7 2 



(o)-t 

P 3 0«1 



(0)~ 
P33°3 



(0)a 

P32 °2 



(26) 



dt 

d 
dt 

d 

dt 
d 

dt 



P3i(t) = g[a 3 pn - p 33 a 3 - p3 2 d 2 + P3oa[] - 7P3i, 



P3o(*) 

P22(*) 
P2l(*) 



r aP30 P(*) + S[«3PlO " P31&1 



(16) 

7P30, (17) 



r a p22 } p(t) + g[a.2Pi2 
g[a 2 pu - P23&3 - P22&2 



P2l4] 



7P22, (18) 



P2o4l 



7P2X) 



(19) 



P21 



P10 



gr a p 



gr a 

7 2 



(o)-t (°)~ 

P20 a l _ P22 a 2 



(0)- 

Poo ai 



(0)-t 

P20 a 2 



(0)- 
P32 «3 



P 3 a 3 



where p\f = pf) 
populations and coherences are real constants 



(27) 



(28) 



py is set based on the fact that the initial 



to express Eq. (11) as 



Now applying Eqs. (26), (27), and (28), it is possible 
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dp 



r a g 2 P33 



2 (0) 

r a g P22 

2 (0) 
r ad POO 
7 2 



2a 3 pa,3 — pd 3 a\ — a 3 d 3 p 



2a\pa,2 — pa2a\ ~ Q-2&\p 



2a\pa\ — pa\ai — a\d\p 



r a 9 2 P32 



2 (0) 



2 (0) 

r a g P20 



2d^pa2 + 2a2/5o,3 — p^a^ — d, 3 d\p — pd.30^2 ~ ^2^3/0 



2aipa 3 + 2a3/5a{ — pa\a\ — a\a\p — a 3 a\p — pa 3 ai 



2a\pa,2 + 2a\pa\ — pa,2ai — a,2a\p — a\a\p — pa\a\ 



(29) 



It is not difficult to observe from the form of this mas- 
ter equation that the two roots of the spontaneous transi- 
tion resemble their counterparts in the three-level cascade 
scheme [3T]. It is good to note that, in addition to the 
correlation pertinent to the cascading process, the initial 
preparation also contributes a new term with a prefactor 

r " 9 1 2 32 since p^ represents the initial coherence associ- 
ated with these atomic energy levels. Based on the fact 
that the cross-correlation terms are indicative of nonclas- 
sical features, it is possible to observe that this master 
equation can be taken as evidence for the existence of 
correlation between the three emitted photons. 

In Eq. (29) there are six different prefactors that are 



not independent altogether, which might make the anal- 
ysis more difficult. As a result, in order to rewrite this 
master equation in a more appealing manner, it appears 
convenient introducing two parameters defined by 



(0) 

P33 I 



(0) (0) 

m = Poo - P22 ■ 



(30) 



(31) 



771 and 772 are basically the population inversions viewed 
from different upper energy levels. In light of the an- 



ticipated initial preparation, 771 and 772 are not entirely 
independent. This can be evinced by the fact that 
p*22 + Poo = 1, which leads to 





P33 



(o) _ 1 + 771+772 

POO — Q ! 



(0) 

P22 



1 + 771 - 2772 



(0) _ 1 + 772 - 2?7i 

P33 — Q 



(32) 



(33) 



(34) 



Moreover, based on the nature of the initial state (Eq. 
pj), that is, ($} = C 3 (0)C 3 *(0), p$ = C 2 (0)C 2 *(0), and 
— Co(0)Cq(0), it may not be difficult to realize that 



TP) 
Poo 

(0) 
P30 



JO) (0) 

P33 Poo ■ 



(0) 

P20 



P22P00 1 ) and P ( 32 



(0) (0) 

P33 P22 ■ 



Furthermore, assuming that the environment modes 
can be represented by a three-mode independent vacuum 
reservoir, it is possible to include its effect following the 
standard approach [43 . In this respect, with the aid of 
the above variable transformation, one readily finds 



dp 
~dt 



= AB 


2a\pa 3 


- pa^al 


- a 3 a\p 


+ AE 


2a^pa,2 


- pa 2 a\ 


~ d 2 a\p - 


\- 2a\pa 3 - 


pa 3 al 


- a 3 alp 


+ AC 


2a\pd,2 


- pa 2 a\ 


- a,2a\p 


-AF 


2aipa 3 


- a 3 dip 


- pa 3 ax -\ 


- 2d\pd\ — 


pa\a\ 


- a\a\p 


+ AD 


2di pa\ 


— pd\d\ 


— a\d\p 


- AG 


2aipa 2 


— pa 2 ai 


— a 2 aip - 


- 2d\pd\ — 


a\a\p 


- P&\ 4 



3 

K x -> 

2 ^ 

i=i 



|^2a 4 /5d- - <i\ ii,fi (Hi)*!, 



iAt At ft 



(35) 



where 



B = 



1 + 772 - 2t?i 
6 



(37) 



2r a g 2 



(36) 



C 



1 + Vi - 2?; 2 
6 



(38) 



G 



D = 



1 + Vi + V2 
6 



(39) 



E _ y/1 - ??i - ??2 + 5?7i77 2 - + rgj ^ 



6 ' 



G = \/l - ??2 + 2?7i - 771772 +vf- 2rj ^ 
6 

and k is the cavity damping constant taken to be the 
same for all modes for convenience. 



IV. DEMONSTRATION OF THE TRIPARTITE 
ENTANGLEMENT 



It is clearly shown that this master equation is essen- 
tially described in terms of A, 771, 772, and k. It is possible 
to infer from the form of the master equation that the 
first two terms to the left indicate the gain in modes 63 
and a.2, in respective order, whereas the third term the 
loss of mode a\. This outcome is fairly consistent with 
the earlier reports on the cascade three-level scheme. As 
already discussed previously, there is a cross-correlation 
between the three modes, primarily in the form of &i and 
0,2, ax and 0.3, and 0,2 and 0,3, separately. Hence from the 
outset, it is not hard to envisage, according to the von 
Loock and Furusawa criteria [TT] , a genuine CV tripartite 
entanglement of a radiation generated by a coherently 
prepared Y-shaped four-level laser. It is obvious that the 
detail of the detectable degree of entanglement depends 
on the strength of E, F, and G, which on the other hand 
heavily rely on A, 771 , and 772 . It goes without saying that 
the strength of the entanglement by and large depends on 
the way the atoms are initially prepared and the rate at 
which the atoms are injected, which believed to give the 
experimenter a considerable freedom for manipulation. 

It is good to note that the relation between 771 and 772 
is so subtle that it can lead to very rich alternatives in 
analyzing the system. For instance, when the atoms are 
initially prepared to be in the lower energy level, one can 
readily see that 771 = 772 = 1. In this case, B = C = 0, 
which indicates that no photon is generated from the two 
upper levels. Moreover, E = F = G = 0, which implies 
that virtually there is no anticipated correlation; as it 
should be. However, assuming the atoms to be prepared 
initially in an equal probability between the three levels, 
Vi =V2 = 0, leads to B = C = D = E = F = G= 1/6. 
This indicates that there is a meaningful correlation 
among the emitted radiations. Basically, these two op- 
tions are the two extreme cases where there is no and 
possible maximum coherence at the beginning, respec- 
tively. 



For the sake of convenience, suppose the atoms are ini- 
tially prepared so that 50% of them to be in the lower 
energy level, that is, 771 + 772 = 0.5 and the remaining 
50% of them are in one of the upper energy levels (let us 
say |3)); 771 = and 772 = 0.5. In this case, one can read- 
ily see that C = F = G = 0, which shows that there is 
only one part of the cascade transitions. Since there is no 
photon with 02, tripartite entanglement is not expected. 
In order to see the situation in depth, with the same 
assumption regarding to the lower energy level, suppose 
the remaining 50% population is equally shared between 
the upper two energy levels, that is, 771 = 772 = 0.25. In 
this case, it is not difficult to assert that all prefactors in 
the master equation are different from zero, which indi- 
cates the possibility for having nonclassical correlations 
among the emitted photons. It is not difficult to observe, 
at this juncture, that a similar outcome could have been 
predicted had the 50% population is arbitrarily shared 
between the upper two energy levels; although the de- 
tails can vary. 

With the same convection, suppose the atoms are ini- 
tially prepared in 50:50 coherent superposition of the up- 
per two energy levels; 771 = 772 = —0.5 (please note that in 
this case the lower energy level is initially unpopulated). 
In this case, one can readily see that D = F = G = 0. 
This indicates that the radiation emitted in the sponta- 
neous transition between the lower two energy levels is 
not correlated with the corresponding upper two transi- 
tions. This can be directly linked to the fact that since 
the upper two energy levels are prepared with a maxi- 
mum coherent superposition between them, the atomic 
transition is basically restricted to transitions from en- 
ergy level 1 3) to |1) and then to \2) and vice versa. The 
chance that the atoms break the established coherence 
and goes over to the lower energy level is quite small. 
That is why even the corresponding mean photon num- 
ber of the radiation which largely depends on a prefactor 
D also can be quite small. In the same manner, it is possi- 
ble to assert that there could be a meaningful correlation 
between the photons emitted during a direct spontaneous 
transition from the upper two energy levels to the lower 
due to the coherence induced by the cascading process 
when the initial coherent superposition is not the max- 
imum possible. The nonclassical correlation in this sys- 
tem seems to show the reminiscent of the absence of the 
bipartite entanglement in a three-level cascade system 
when the lower and the upper energy levels are prepared 
in a maximum coherent superposition |32j . 

One may deduce from this interpretation that the non- 
classical correlation that leads to a tripartite entangle- 
ment can be induced by initially preparing atoms coher- 
ently between the upper two energy levels and the lower 
with arbitrary, other than zero, probability. It is also 
envisioned that the cascading transitions are a very vi- 
tal mechanism in connecting the upper two energy levels 
with the lower. Therefore, no doubt that, properly har- 
nessing the utility accorded with the initial preparation 
and cascading mechanisms results a genuine tripartite 
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entangled light. Based on earlier studies in the three- 
level cascade laser [361 EI], it is equally expected that 
externally pumping the atoms can establish a coupling 
between energy levels in which direct spontaneous tran- 
sitions are dipolc forbidden that can significantly improve 
the tripartite nonclassical correlations. 



V. CONCLUSION 

The detailed derivation of the master equation that de- 
scribes the radiation emitted from a coherently prepared 
nondegenerate ^-shaped four-level correlated emission 
laser is presented. In view of limiting the otherwise aris- 
ing complications, the atoms are presumed to be pre- 
pared in an arbitrary perfect coherent superposition of 
the three energy levels, where the intermediate energy 
level is taken to be unpopulated at the beginning to make 
use of the adiabatic elimination technique. In setting up 
the laser, the initially prepared atoms are presumed to be 
injected into a triply resonant cavity. To pave the way for 
further in depth analysis, the approaches by which the 
corresponding stochastic differential equations can be ob- 
tained and the resulting equations are solved are outlined 
in the Appendix. Moreover, the rate equations that de- 
scribe the time evolution of various correlations that can 
be required in the study of the quantum features and 
statistical properties of the radiation and the procedure 
applied in solving them are provided. 

It turns out that the quantum system under consider- 
ation can be a source of a continuous variable tripartite 
entangled light under certain conditions. Detailed inves- 
tigation shows that due to the cascade transitions, the 
emission-absorption mechanism which is guided by the 
induced coherent superposition is found to be responsi- 
ble for establishing the required correlation between the 
emitted photons. Further analysis based on varying the 
way the atoms are initially prepared shows that coher- 
ently coupling the three atomic energy levels is very vital 
in generating a tripartite entangled light. It is unequivo- 
cally asserted that leaving one of the upper energy levels 
unpopulated at the beginning leads to a bipartite entan- 
glement at best since there is no way for the atoms to 
go to the unpopulated level in the course of the process. 
However, initially preparing the atoms in the upper two 
energy levels other than in a possible maximum coher- 
ence leaving the lower energy level unpopulated can lead 
to the appearance of the tripartite entanglement since 
the upper two energy levels can be coupled to the lower 
via the cascade transitions. 

In relation to the similarity of the result of preparing 
the atoms in the coherent superposition of the upper two 
energy levels with the corresponding three-level scheme, 
it is expected that an external driving mechanism can be 
used to improve the generated entanglement in some re- 
spect. Furthermore, it may not be hard to realize that a 
highly intense light can be generated since the injection 
mechanism allows to send as many as required atoms 



through the cavity over a longer period of time without 
significantly exposing them to fluctuations and broad- 
ening associated with heating. Hence, this study by and 
large tries to show that the nondegenerate ^-shaped four- 
level scheme can be a source of reliable bright genuine tri- 
partite entanglement with much more promise of flexible 
arrangement essentially by combining the initial prepara- 
tion and external driving options. Even though starting 
with the form of the master equation and the values of 
the involved prefactors yield encouraging outcomes, it is 
incontestable that in depth analysis is still lacking. Ow- 
ing to the involved rigor, length of the expressions, and 
complications of the different viable scenarios, a simpler 
and more specialized analysis is deferred to subsequent 
communications . 
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APPENDIX 

In order to study the quantum features and statistical 
properties of the radiation, different auto-correlations, 
that refer to the pertinent mean photon numbers, and 
cross-correlations are required. In many instances, it is 
possible to obtain these correlations by making use of 
the master equation and then solve the resulting coupled 
differential equations. To this effect, employing Eq. (35) 

and the fact that f t {6{t)) = Tr (j£d) , where 6 is an 



operator, the following essential correlations are obtained 
(Please note that the full list is unnecessarily long) 



dt 



a\(t)) 



AD\(a\(t)) 



dt 



+AF{a 3 (t))+AG{a 2 (t)), 
~{a 2 (t)) = -(*-Ac) (a 2 (t)) 

+AE(a 3 (t))-AG(a\(t)), 
j t {a,{t))=-(^-AB) (a 3 (t)> 

+AE(a 2 (t))-AF(a\(t)}, 
(ala 3 ) = (2AB - n)(dla 3 ) + AE[{a\a 2 ) + (a 3 4)] 
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dt 



dt 



-AF[(ala\) + (a 3 oi)] + 2AB, 



(a\a 2 ) = (2AC - n)(d\a 2 ) + AE[(a\a 2 ) + (a 3 a\)} 



-AG[(a\a\) + (a 2 a 1 )}+2AC\ 



(ajai) = -(2 AD + ^(ajai) + AF[(ala\) + (a 3 fli>] 



+AG[(aIa 2 > + (a 2 ai)] - 2 AD, 



where 




—AG 




-AC 


—AE 


-AE 


l-AB 


1 _ * 






I 

1 ' 


v «3(t); 








/ 2 (*) 


1 


\h(t)j 





(AA) 



— (d\a 2 ) = (AB + AC- K){a\a 2 ) - AG{a\a\) 



+AE[(a\a 3 ) + (a\a 2 )] - AF(a 2ai ) + 2AE, 



dt 



(a 3 ai) = (AB - AD - /c)(a 3 ai) + AG{a z a 2 ) 



+AF[(a%) - (al^)] + AF, 



dt 



(a 2 oi) = (AC - AD - n)(d 2 ai) + AE(dia 3 ) 



+AF(a 2 a\) + AG[(a\a 2 ) - (&%)] + AG. (Al) 

Assuming that c-number expressions associated with 
normal ordering are mathematically more appealing, the 
first three equations are rewritten as 

±a* 1 (t) = -(*+AD) al(t)+AFa 3 (t) 



+AGa 2 (t) + f*(t), 



dt 



a 2 (t) = - (| - AC) a 2 (t) + AEa 3 (t) 



-AGa\(t) + f 2 (t), 



dt 



o 3 (t) = - (f - «3(*) + AEa 2 (t) 



-AFa\(t) + h(t), 



(A2) 



where fi(t)'s are the pertinent stochastic noise forces. 

For the sake of convenience, these equations can be put 
in a more compact from as 



dt 



TZ(t) = -MTZ(t) +Af(t), 



(A3) 



In principle, this coupled differential equations can be 
solved following a somewhat lengthy but straightforward 
algebra. First of all, it is desirable and possible to diag- 
onalize the matrix A4 using the eigenvalue equations in 
which MVi = XiVi, where Ws are the eigenvectors and 
Ai's are the corresponding eigenvalues. For a 3X3 matrix, 
although the involved rigor is lengthy, it is possible to find 
both the eigenvalues and the corresponding eigenvectors, 
with the property that VV" 1 = 2 and V = V^MV, 
where V -1 is the inverse of the matrix constructed from 
relevant eigenvalues and V is the diagonal matrix corre- 
sponding to Ai. With this arrangement, the solution of 
Eq. (A3) can be proposed as 

■R(t)=[Ve~ vt V~ 1 ]ll({))+ I [Ve- v{t - t ' ) V- 1 ]N(t')dt l . 

Jo 

(A5) 

It is common knowledge that 72.(0) describes the prop- 
erties at the beginning of the lasing process. If initially 
the cavity is assumed to be in a three-mode vacuum state, 
one can readily disregard the contribution of the first 
term in Eq. (A5) and it is also possible to verify by tak- 
ing the expectation values of the expressions in Eq. (A2) 
and comparing them with the first three expressions in 
Eq. (Al) that 



(A/-(t)> = 0, 



(A6) 



which implies that (fi(t)) = for i = 1,2,3. This essen- 
tially reflects the stochastic nature of the noise. 

Furthermore, applying the various rate equations, dis- 
tinct terms of Eq. (A5), and assuming the noise force at 
later time does not affect system variables at earlier time 
result in 



F(t',t") 




t"), (Al) 



where T(t,t') = (Af(t)Af T (t')), in which T stands for 
complex conjugate transpose. Based on the results ob- 
tained so far, it can be observed that 



(K(t)K T (f)) 



o Jo 



T(t,t')F(t',t")T T (t,t")dt'dt'' 



(A8) 
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where V(t,U) = Ve^*"**) V" 1 . 

In principle, carrying out the involved matrix manip- 
ulations and then term by term integrations yield the 
correlations required for studying the quantum features 
and statistical properties of the radiation without further 



approximation. Even though the number of terms to be 
handled is somewhat large, carrying out the integration 
is over simplified by the presence of the (5-function asso- 
ciated with the correlation of the noise forces (F{f ,t")) 
and the exponential dependence. 
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